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J> . Abstract: We study the factorization and resummation of s-channel single top quark produc- 

•rH , 

tion in the Standard Model at both the Tevatron and the LHC. We show that the production 
cross section in the threshold limit can be factorized into a convolution of hard function, soft 
function and jet function via soft-collinear effective theory (SCET), and resummation can 
be performed using renormalization group equation in the momentum space resummation 
formalism. We find that the resummation effects significantly reduce the factorization scale 
dependence of the total cross section at the Tevatron, while at the LHC we find that the 
factorization scale dependence has not been improved, compared with the NLO results. 
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1. INTRODUCTION 

Recently, the analysis of data by DO jjj and CDF collaboration has confirmed the obser- 
vation of single top production at the Tevatron. Due to the fact that its expected production 
cross section is small ( in combined s and t channel, <r s t ~ 2.9 pb ||) and the processes lack 
significant signal feature, such discovery signifies a success in both experiment and theory, 
and provides good opportunities for testing the Standard Model (SM) and searching for new 
physics. 

As we know, at hadron collider, single top quarks are produced via three different modes: 
s-channel (qq' — > bt), t-channel (bq — > tq' and bq — > tq'), and associated tW production 
(bg — > tW~), each sensitive to quite different physics H|. First of all, since all three production 
modes are directly proportional to the CKM matrix element | V^| 2 , a measurement of the cross 
section provides a unique direct probe to Vtb, and can constrain models with fourth generation 
quark. The s-channel single top production is rather sensitive to the interaction mediated 
by extra heavy particle and the uncertainty of partonic luminosity for this mode is relatively 
small. Thus although this mode has the least cross section, it's a very important channel 
in searching for new physics. The t-channel is sensitive to physics which modifies top decay 
properties, and the associated production channel can be a good measurement of the W — t — b 
vertex. Thus a precise understanding of production cross sections and theoretical uncertainty 
are important. In particular, higher order QCD corrections to the cross section are necessary 
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to improve theoretical predictions. There have been a lot of NLO calculations of the single 
top production in the literatures [§, g, 0, |, |, |, 0, 0, 0, 0, ||, H, |§ , with or without 
subsequent top quark decay. Recently, implementations of these results into NLO Shower 
Monte Carlo (MC@NLO or POWHEG) also appeared ||, g]§ 



Furthermore, the consideration of phase space logarithms in higher order QCD effects, 
which origin from incomplete cancelation of real soft gluon emission and virtual corrections, 
are also important to the theoretical predictions. They often dominate the hadronic cross 
sections and can be systematically resummed to all orders in perturbation theory j pO| , pl{ . For 
the single top production they have been calculated to NNLO at Next-to-next-to-Leading- 
Logarithmic (NNLL) accuracy |2^, ^3|, but with some important NNLL logarithms omitted. 

In recent years it has become popular to use SCET to resum the large phase space 
logarithms. SCET is an effective field theory describing the QCD interaction between collinear 
and soft particles [ 24 , 2^. 26 , ^?j] , which can correctly reproduce the long distance behavior of 
QCD, while the short distance information is encoded in the Wilson coefficients from matching 
the full theory to SCET. In the past decade, SCET has proved its usefulness in describing high 
energy hard scattering processes. These include deep-inelastic scattering 
1 34, 35], Higgs production [36, 33, |37 



,30,81, Drell- 



39 



e + e annihilation to 



Yan production |3^, 
hadrons f4|, H> E§ E 



color-octet scalar production |Q, direct photon production [47], 
direct top quark production via FCNC coupling and top quark pair production g|, HI H HI • 
In this paper, we will further study the threshold resummation effects on the s-channel 



production cross section to all orders in QCD at NNLL accuracy in momentum space ]30j, 
utilizing SCET. First of all, we show that the total cross section for s-channel single top 
production can be factorized schematically as 



a 



f®f®H®S®J, 



1.1) 



where / is the initial state nonperturbative parton distribution functions (PDFs); H is the 
hard function, which encodes the short distance interaction information; S is the soft function, 
which describes the soft correlation between different color objects; and J is the jet function, 
describing the final state collinear emission associated with the jet. The factorized cross 
section we derived is valid in the hadronic threshold. In this limit, the partons initiated the 
hard scattering carry almost all of the hadron momentum, and the final state configuration 
consists of a top quark, a narrow hard jet and the remaining soft radiations. In this work, we 
are only intereseted in the inclusive total cross section, thus do not consider the top quark 
decay effects. For this reason, our results can not apply to the isolated s-channel single top 
cross section measured at the Tevatron, where explicit experimental cuts on the final state 
leptons and jets are required, but only serve as part of the total single top production cross 
section. When combined with the results of t-channel and tW associated production channel, 
our results can provide the most accurate perturbative predictions for total cross section of 
single top production. 

This paper is organized as follows. In Section |2[ we briefly review the basic ingredients of 
SCET. Section |3| discusses the kinematics at threshold. In section |]we derive a factorization 
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formula for the resummed cross section in momentum space. We give a NNLO expansion of 
our resummed cross section in section ||. Section |(| contains a brief numerical discussion and 
we conclude at section |7[ 

2. Brief introduction to SCET 

To describes collinear field in SCET, it is convenient to define a lightlike vector = 
(l,n),n 2 = 1. Any four-vector can be decomposed with respect to and n M = (1,— n) 
as 

^ = ry+! + y+C, (2.1) 

with l + = n-l and l~ = n-l. The momentum of a collinear particle moving along the n^ 
direction has the following scaling 

I? = (p + ,p-, Pnl .)~(\ 2 ,l,\), (2.2) 

while for a soft particle, the momentum scales as 

g~(A 2 ,A 2 ,A 2 ), (2.3) 

where A <C 1 is a small expansion parameter in SCET. E.g., for an energetic jet with invariant 
mass mj and energy Ej, A = mj/Ej. From the momentum scaling, one can see that the 
interaction between collinear fields of different directions rii and nj with m ■ rij 3> A 2 will 
inevitably change the momentum scaling, thus is forbidden in SCET, but can be included as 
an external current in our computation. The soft fields, on the other hand, can interact with 
any collinear field without changing the scaling. 

A n-collinear quark and gluon field can be written as 

Xn{x) = Wl{x)£ n (x), 

(2.4) 



W^ ± W n (x) 



where 

is the collinear covariant derivative and the label operator V is defined to project out the 
large momentum component of the collinear field, e.g., VnCn = P^£,n- Here we have split p 
into a sum of large label momentum and small residue momentum, 

p»=p» + k>*, with =P~Y +P n±- ( 2 - 6 ) 

The collinear Wilson line, 

W n (x) = Pexp (ig a J dsn-A a n {x + sn)t a ^ , (2.7) 
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which describes the emission of arbitrary n-collinear gluons from an n-collinear quark or 
gluon, is constructed to make the collinear fields as defined in Eq. ( |2,4D invariant under the 
collinear gauge transformation. The operator P is the path-ordered operator acting on the 
color generator t a . 

At the leading order in A, only the n-A s component of soft gluons can interact with the 
n-collinear field. Such interaction is Eikonal and can be removed by a field redefinition [27]: 

Xn (x) = Y n (x) X ^(x), 

B^{x) = Y n (x)B^(x)Y^x), (2.8) 

where 

Y n (x) = Pexp (ws J dsn-A a s {x + sn)t a ^j (2.9) 
for an incoming Wilson line [p2|| . And for an out going Wilson line, it is defined as 

Y n (x) = Pexp (-ig s dsn- A a s {x + sn)t°^j . (2.10) 

The fields with superscript (0) now are decoupled with soft gluons. Without confusion, we will 
neglect the superscript below. After the field redefinition, the leading order SCET lagrangian 
is factorized into a sum of different collinear sectors and soft sector, which do not interact 
with each other, 

£ S CET = E £ n? +£, + ••• • (2-11) 



3. Analysis of kinematics 

In this section, we introduce the relevant kinematical variables needed in our analysis. As 
an example, we consider the subprocess u + d — > t + X. The subprocesses induced by gluon 
splitting are power suppressed in the threshold limit, therefore will not be considered in this 
paper. First of all, we define two lightlike vectors along the beam directions, n a and rif,, which 
are related by n a = fib- Then we introduce initial collinear fields along n a and to describe 
the collinear particles in the beam directions. In the center-of-mass frame of the hadronic 
collision, the momentum of the incoming hadrons can be written as 

^ = £ C My, If = Ea A -±. (3.1) 

Here Eqm is the center-of-mass energy of the collider and we have neglected the masses of 
the hadrons. The momentum of the incoming partons, with a light-cone momentum fraction 
of the hadronic momentum, are 

p a = x a E C M-^-, Pb = XbE C M-Jr- (3.2) 
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At the hadronic and partonic level, momentum conservation means 



P a + Pb = q + Px, (3.3) 

and 

Pa+Pb = Q + Px, (3.4) 

respectively, where q is the momenta of the top quark. We define the partonic jet with jet 
momentum px to be the set of all final state partons except the top quark in the partonic 
processes, while the hadronic jet with jet momentum Px contains all the hadrons as well as 
the beam remnants in the final state, except the top quark. Explicitly, px = pi + k, where p\ 
is the momentum of the final state collinear partons forming the jet and k is the momentum 
of the soft radiations. Such division of momentum is artificial and we have to integrate over 
the soft momentum to obtain a physical observable. We also define the Mandelstam variables 
as 

s = {P a + P b ) 2 , u={P a -q)\ t ={P h -q) 2 (3.5) 



for hadrons, and 



s = (p a +Pb) 2 , u=(p a -q) 2 , t=(p b -q) 2 (3.6) 



for partons, respectively. In terms of the Mandelstam variables, the hadronic and partonic 
threshold variables are defined as S4 = P\ = s + t + u — m 2 and S4 = p\ = s + t + u — m 2 , 



where rat is the mass of top quark. The hadronic threshold limit is defined as S4 — > [53]. 
In this limit, the final state radiations and beam remnants are highly suppressed, leads to a 
configuration consists of a top quark and a narrow jet, as well as the remaining soft radiations. 
Taking this limit requires x a — > 1, Xb — > 1, S4 — > simultaneously. Thus, the hadronic 
threshold enforces the partonic threshold. However, the reverse is not true. The partonic 
threshold 54—7-0 does not forbid a significant amount of beam remnants. We note that in 
both hadronic threshold limit and partonic threshold limit, the top quark is not forced to be 
produced at rest, i.e. it can have a large momentum. For later convenience, we can also write 
the threshold variable as 

Si=Px = (Pa+Pb ~ q? =p\ + 2k+E 1 + 0(k 2 ), (3.7) 
where k + = nyk, E\ is the energy of the quark jet and n\ is the lightlike vector associated with 



the jet direction. Note that our definition of s\ is different from [22|, in which the definition 
S4 = {p a + pb — pi) 2 — m 2 is adopted, where we put a bar on 54 to distinguish the definition 
in |z| from ours. We point out that the meaning of such choice is not clear, because S4 
doesn't vanish when there are collinear gluon emitting from the final state b-quark. However, 
as we know, there are large logarithms associated with such collinear gluon emission, thus, 
the definition adopted in p^] could miss some large contributions. 



In the threshold limit (54 —> 0), incomplete cancelation between real and virtual correc- 
tions leads to singular distributions a"[ln m (s4/m^)/s4]+, with m < In — 1. It is the purpose 
of threshold resummation to sum up these contributions to all orders in perturbation theory. 
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4. Factorization in SCET 



In this section, we derive the factorized cross section formula for s-channel single top produc- 
tion in SCET, following the convention and formalism of |54|, [55|]. With appropriate changes, 
our formula can be extended to the resummation of t-channel and tW associated production 
channel as well, which will be presented elsewhere [?, |57p . 

Total cross section for s-channel single top production can be written as [54]: 

_. res. 

<r= 5 i5-X;<I|Oi(0)|X>(X|O ie (0)|i)(27r) 4 d 4 (P o +n-px) 

Z£j CM x 

= / d 4 x(I\OUx)\X)(X\O x (0)\I) 

Z£j CM x J 

1 f , 4 f d A k ik . x f d 4 p 



9 f?2 

Z - C/ CM x 

1 res - r I I 

1 ^-^ / d p 



res. „ ,4 

E j j^(I\OH0)\X)(X\O(p)\I), (4.1) 



2E 2 

zr/ CM x 

where |J) = \P a Pb) denotes the initial state protons (anti-protons), O x (x) is the relevant 
operator responsible for the underlying hard interaction, and 0(p) its fourier transformation. 
Here we distinguish the position space operator from the momentum space one by a subscript 
x. The restriction on the sum over final states \X) is that we include final state configuration 
consists only of a top quark jet whose 3- momentum is in the direction of ni, an anti b quark 
quark jet in the direction of ni, and soft radiations. This is the configuration that is relevant 
to threshold resummation and we are interested in. Under this condition the final state can 
be written as \X) = \XtXiX s ) , where \Xt), \X\) and \X S ) denote top quark jet, anti b quark 



jet and remaining soft radiations, respectively. In the second line of Eq. (4.1), we have used 



the momentum conservation delta function to shift the operator O x to point x, and in the 
third line written the operators in momentum space, which are then matched onto SCET 
operators. We chose to match directly in momentum space, which has the advantage that the 
cumbersome label summations and residual soft integration can be combined to a full four 
momentum integral |)4|] . After matching we obtain 

i \ f d 4 p a d 4 p b d A pi d A p 2 d 4 k s 
° {P) = J (2tt)M2^ (2vr)4 (2^ frff ^™^*) 

x O in ( Pa , Pb )0 0Ut (pi , p 2 ) O s ,i (k s )(27r) A 5^(p-p a -p b + Pl +p2 + k s ). (4.2) 

In the above equation, we have made all the spin, Lorentz and color indices implicit. For 
example, at the LO, the hard matching coefficient Ci reads 

C i = i ■ 2n T Kd \ ^ (7 M (1 - 7 5 )) 75 (7„U " 7 5 )) a/? 5n, (4.3) 



2 sin 2 0w((Pa + Pb) 2 — M } 



2 

XV > 
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where a is the fine-structure constant, Vij the CKM matrix element, 9w the electroweak 
mixing angle, and M\y the mass of the W boson. Here we have chosen the s-channel singlet- 
octet basis as the independent color structure for this process: 



|ci) = 5 cd S e f, 



C2 



(t a )cd(t a ) e f, 



(4.4) 



and I = 1 or 2 is an index in this color space. Thus the kronecker delta function in Eq. 
can be understood that at the LO the final state tb pair of s-channel single top production is 
a color-singlet. In Eq. fl4.2| ) Oi n denotes the effective operator responsible for annihilating an 
initial state collinear up-type quark with momentum p a and an down-type anti-quark with 
momentum p^, which can be written as 



0%),in = X C a(-Pb)X$(.Pa), 



(4.5) 



and 0U t is the effective operator responsible for the creation of final state anti-bottom quark 
with momentum pi and top quark with momentum p2 , which can be expressed as 



O 



fi 



(4.6) 



Note that we have taken the mt — > oo limit at fixed top jet radius and then described the 
top quark in terms of the heavy quark effective field [58] with label velocity v. The fields 
in Eqs. ( |4.5| ) and (|4.6j ) are defined with field redefinition in Eq. ( |2.8[) , thus they no longer 
interact with soft degrees of freedom. 

The soft operators Os,i, which are responsible for the soft interactions between different 
collinear sectors and top quark, are expressed as 



OsAk s ) 
OsAks) 



d xe 



T 



d 4 xe~ ik °- x T 



yl(x)Y na (x)) C (Yl(x)Y ni (x)Y 
Yl (x)t a Y na (X)) ^ (?i (x)t a Y ni (x) 



(4.7) 



where the time-ordering operator T is required to ensure the proper ordering of soft gluon 
fields in the soft Wilson line. 

Using the notation $2 = {Pa,Pb',Pi,P2} to express a phase space point with d&2 = 
d A p a d A pbd 4: pid 4: p2/ (2tt) 16 and $2 — k s = p a + Pb — Pi — P2 — k s , we can write Eq. (4J) in a 
compact form 



CM x 



x{I\{0[ n O' mt O' s J\X t X 1 X s ){X t X 1 X s \{O in O wt O SiI )\I). 



(4.8) 



As we mentioned before, different collinear sectors are decoupled due to field redefinition, 
and thus the matrix element in Eq. (|4.8j ) can be factorized into a product of several matrix 
elements, which obey certain renormalization group (RG) equation. 
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In the following, we further show the matrix elements mentioned above. First, we deal 
with the top quark sector. Since we have decoupled the soft interaction by field redefinition, 
the top quark now should be regarded as a non-interacting particle, which can be written as 

£ / 1% (p| ^ {p ' 2) ]xt) {xt 1 ^ fe) |o) - 



x 



where summation over final state \Xt) gives rise to a top quark phase space integral. Next, 
we define the soft function by the soft matrix element as 

d k e d kc 



S{k+ - n v h,]\X,){X,\0^/{k,)\<S), (4.10) 
where we have inserted into the above equation an identity operator 

dk + 5[k + -n r k s }. (4.11) 



Note that the summation over final state can be performed Ylx s l^s)(^s| = 1 since there is 
no restriction in the summation when written between the soft operator and also there is no 
explicit dependence on \X S ). Since we are only interested in the total cross sections, the final 
state top quark, jet function and PDFs can be considered to be diagonal in color space. Then 
we can contract them to obtain the soft function matrix 

Sjj(k + , y) = 5f , fS e ' e 5 d ' d 5 c ' c S f J ' I e ' d ' c ' cdef (k + , y). (4.12) 
At the LO, it can be written as 

( N 2 \ 

S(fc+, M ) = *(*+) ( Q c j • (4.13) 

At the NLO, the calculation of soft function boils down to the evaluation of eikonal dia- 
grams p^] . Since the virtual corrections in SCET vanish, only real emission diagrams, which 
are shown in Figs. |l|, are needed to be evaluated. The detail calculation of these diagrams 
are given in the appendix. 



For the final state anti-quark jet sector, we have [54] 



= sn lm>^)j (AM) - (414) 
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(a) (b) (c) 




(d) (e) (f) 




(g) 

Figure 1: Diagrams contribute to the soft function, the double line represents the top quark line. 
Diagrams (a), (c), (d), (e), (f) and (g) have their mirror image diagrams. There are also diagrams 
with soft gluon connecting to the same massless quark line on the two side of the cut, which are not 
shown here because they vanish identically. The cross dot in the diagrams represents either singlet or 
octet operator. 

Again, summation over collinear state has been performed and J is the spin and color singlet 
jet function, which can be defined as 

e(p°)p-J( P 2 ) = -^r/ 0^Tr(O\x(- P 'nix(-p)\O). (4.15) 

Finally, the initial state collinear sector reduces to the conventional PDFs p^ ], of which the 
matrix element for n a direction is 

/ |£||<^to4Mip.> = ±f<[ f Xcm^ /(*„,, X4,6, 

and similarly for the matrix element for n& direction. Thus the momentum of incoming 
partons are given by p a) b = x a ,bEcMn^ h /2. 
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Combining the above expressions, we obtain (up to power corrections) 
1 If 1 dx a dxb f d 3 q 



a 



2££ M 4iV c 2 Jo x a x b 
x J dk + Sj^kf , /i)(2vr)J( S4 - 2k+E u fj,) 



f d A q 



with 



A, 



0,ij 



\Vij\ 2 \V tb \ 2 



(i-m 2 t )t 



(4.17) 



(4.18) 



am* iv-" {s-M^r 

At the LO, the hard function Hjj is normalized to 5n6ji. In general, it is related to the 
amplitudes of full theory by (5^] 

1 



A 0,ij H IJ 



(ci\ci)(cj\cj) 

1 

{ci\ci)(cj\cj) 



(ci\M<®)(MW\cj 



(0)| 



(dlMgDiMilllcj) + (cj\M^)(Mgl\cj)) , (4.19) 



where |-A4 ren ) are obtained by subtracting the IR divergences in the MS scheme from the UV 
renormalized amplitudes of full theory. To 1-loop order, it reduces to evaluating in the full 
theory the 1-loop Feynman diagrams in Fig. ||. The complete 1-loop hard function is shown 
in the appendix. 






+ c.t. 

Figure 2: 1-loop Feynman diagrams for the hard function. The double line represents the top quark. 
The hard function is a 2 x 2 matrix in the color space. The RG equation it obeys reads 



din fj, 



r^H + Hr^. 



(4.20) 



The relevant anomalous dimension matrix is process dependent and can be expanded in 
|k, with the relevant expansion coefficient given in the appendix. The explicit form of the 
anomalous dimension matrix for s-channel single top production can be extracted from a 
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more general result given in the Ref. |]59[| . Explicitly, for the independent color basis we have 
chosen, it is given by 



rax 



77CF7cus P (a s ) hi — j + lh{a 



1 + T h , 



with 



and 



, v „ / \ , s(s m?) 
jh{a s ) = G F j cusp (a s )ln — A h 7ft («s), 



r/i = 7cusp(a s ) 







C F i u{u-m? t ) 
2N C 111 t(t-m 2 t ) 



, u{u-m4) _/Vc i u(u-mj) J, u{u-mj) 

111 t(t-m?) 2 m s(s-mf) £(i-m t 2 ) 



(4.21) 



(4.22) 



(4.23) 



The expression for 7^ and 7 CU sp can be found in the appendix. Note that we have only 
retained the real part of the anomalous dimension matrix. We have checked that the 1-loop 
hard function exactly obeys this equation, as it must. Details of it are presented in the 
appendix. 



The solution of the RG equation shown in Eq. ( 4.20 ) can be obtained by diagonalizing 



the anomalous dimension matrix [60|. A detail example for how to do this in effective 
field theory can be found in the Ref. pi]]. To NNLL accuracy, the result can be expressed as 



irr, 



Hr(/x) = exp [6S D l(m/i>m) ~ 2a h(Hh,n)} \ 



(4.24) 



with 



where 



Ur(^,^) = (1 + 



4ir 



7 - h r° 



The Sudakov exponent are given by 

SclC^aO = - 

ar(u,ij) = - 
a h {v,n) = - 



2/? 2 " 2/3 



asM dA CF7cus P (A) 

^ dA CF7cUSp(A) 



/3(A) 



a s (V) 
a s (u) PW 



47T 



J), 



(4.25) 



(4.26) 



(4.27) 
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Next, we discuss the jet function. Unlike the hard function which obeys a local RG 
equation, the jet function J(p 2 ,fi) satisfies a RG equation which is non-local in p 2 [61]: 



dhx/j, 



P 2 \ 
A 4 / 



+2^ 7cusp (a s ) f d g 2 J(P '^" J 2 (g ' /i) . (4.28) 
Jo P ~ 1 

This equation is solved with the help of Laplace transformed jet function : 

~ O 2 f°° v 2 

j(\n^^)= dp 2 exp(--f— )J(p 2 ,fi), (4.29) 
M Jo Q z efE 

which satisfies the RG equation 



(iln/i 



O 2 \ ~ O 2 

-2C F7 cus P (a s ) ln^- - 2 7i (a s )J j(ln \,//), (4.30) 



and can be easily solved now. The solution, after transformed back to momentum space, 
is @ 

J(p ,fi) = exp(-4S , DL(Aij,Ai) + 2a j (nj,ti))j(d Vj ,[j l j)-jp I ^ I ( 4 - 31 ) 

where i]j = 2ar(/^/,/i). 

Finally, we need the RG equation of the soft function, which can be obtained by noticing 
the fact that the hadronic cross section in the threshold region should be independent of the 
arbitrary scale [i. Schematically, 

^-/0/®H«S0J = O. (4.32) 
In [i 

Based on this fact, we have 

s(]n^,^=^ dfc + exp(--J^) S(fc+,M), (4-33) 

with 



where 



T s = (c F j cusp (a s ) In ^ - %(a s ij 1 + (4.35) 



j s (a s ) = -C F j cusp (a s ) In 2n » n&?; ni + 7s ( as ) ) (4.36) 

n a -n\ n b -ni 
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and 7 S is given in the appendix. Note that in obtaining the RG equation Eq. (f4.34j) , we have 
used the DGLAP evolution for the PDFs in the x — > 1 limit: 

df q/N {x,fi) f 1 dz f q/N (x/z,fi) 

— ^— - — = 2j 4> {a s )f q/N (x,n) + 2C F 7 cusp (a s ) J - _ . (4.37) 

We stress that Eq. ( [4.34j ) is derived entirely from RG invariance of the resummed cross section 
in the threshold limit. If the explicit 1-loop soft function do obey this equation, it serves as 
a non-trivial check on the RG invariance of our result. We show in the appendix that this is 
indeed the expected. The solution of the RG equation in momentum space is 

S R (k + ,fi) = exp[-2S DL (n s ,n) - 2a s (fj, s , fj,)] 

+ _ 1 fk + \ Vs e~ lEris 

,V R (,,, s )s(d Vs ,» s )lJ R (»,» s )-l y , (4.38) 

where rj s = 2ar(^ S iA i )- Combining the above ingredients, we obtain the resummed cross 
section for s-channel single top production 

thres V- ^ f 1 d Xa dx b f d?q 

x / dk + Tv[H R {fi)S R (k + ,/j)}J( S4 -2E 1 k + ,fi), (4.39) 

J o 

where we have included a summation over different partonic channels. 



5. NNLO expansion of resummed cross section 

In the traditional approach to threshold resummation, the evolution equations for the factor- 
ized cross section are solved in Mellin moment space rather than momentum space. It has 
been demonstrated that in Drell-Yan production, the two approach are equivalent up to l/N 



corrections when making the scale choices jj,h = M and /j, s = M/N [34]. Here N is the Mellin 
moment and M is the mass of the Drell-Yan lepton pair. Note that when N is very large, 
the running coupling constant a s (fj, s ) blows up. Furthermore, in order to obtain the physical 
cross section in momentum space, one needs to invert the Mellin transformation numerically, 
which is ambiguous since the expression to be inverted has a Landau-pole for large N, al- 



though there have been several prescriptions for dealing with the Mellin inversion |6^, 63 , |64j . 
On the other hand, in the momentum space approach, we do not encounter such problem 
because fj, s is constrained to be well above the Landau-pole singularity. 

It has also been advocated that the resummed cross section in Mellin moment space 
can be used as a generator of the fixed order perturbation expansion [^5|. The ambiguity 
due to different Mellin inversion prescriptions can be avoided because at the fixed order no 
prescription is needed to invert the Mellin space results. In this way, partial NNNLO threshold 
singular terms were obtained for s-channel single top production |36|, ^37], [22fl , 
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In order to obtain the fixed order expansion in momentum space resummation formalism, 
we rewrite the cross section in terms of integral over singular distributions of S4. In the center- 
of-mass frame of the top quark and the recoiling jet, we can parametrize the momentum of 
p a , pb and q as @ 



Pa = —(1,0,0,1), p b = —(1,0,0,-1), q= (E q ,0, |q|sinx, |q|cosx), (5.1) 



where 



E„ 



i+u — 2m? 



u — t 



-, cosx 



2 ^ ^{i + u- 2m 2 ) 2 - Asm 2 

Hence the phase space measure of the top quark is 

d 3 q 



(5.2) 



2£g(2vr) 3 8vr 2 



dE q dx s'mx\ E 2 - m\ 



1 



167T 2 S 

and the cross section can now be rewritten as 



-dtdu, 



Jm^/s Jt x a Jm^—s JO 

where sf 18 "* = s + t + m 2 s/(t — m 2 ), and 
da thres 



dSA fi/N a { X ^ VF)fj/N b {T/x a , Hf) 



da 



thres 



dtdu 



(5.3) 



(5.4) 



dtdu 



U-kN 2 s 2 J 



dk + Tr[H i? (/x)S R (fe+, fi)] J(s 4 - 2E 1 k + , f x). (5.5) 



Following [34 1, we derive the threshold singular distributions by setting Hh, fJ, s and fij equal 
to the common scale fj,, which is conveniently chosen as the factorization scale [ip. In the 
following, we show all the threshold singular distributions up to NNLO, 



64irs 2 d 2 a^ nd a s 

-= 6(s 4 ) + ^ (A 2 D 2 + A X D X + A 5(s 4 )) 



Xq ijdi du 



4tt 



^ j (S4D4 + B 3 D 3 + B 2 D 2 + B\D\ + B 5(s 4 )) 



where 



Si 



mi 



(5.6) 



(5.7) 



is the conventional plus distribution, and its integral with a regular function /(S4) is defined 
as 



ds A f(s A )D n = I ds 4 [f(s 4 )-f(0)]-\n n - 1 ^. 
Jo Si rnf 



(5.8) 
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The coefficients A n and B n are given by 



C Fl ° cnsv {L h + 2L s ) + 1 j -2f s , 



g| + 4 + 4 + C F7c ° usp (-^ + L» - '- ) + ^(7° 



A 2 = 3Cf7cus P > 

B 3 
B 2 



T 2 TT 2 



7^) 



27% 



§cf( 7 ° 



2 

CUSp y ! 



cusp 



(9Cf7c°us P ^^ + ^C Fl ° cusp L s + 9 7 ° - 18 7s ° - 5ft) , 



CW° 



cusp 



3c 



(5.9) 
(5.10) 

(5.11) 

(5.12) 
(5.13) 



ii 



r 2 



+ 3c h u + 3cj - L h (3j° h - 5 7 ° + 4 7 ° + /3 ) + 4 7 °L S - 14 7s °L s (5.14) 



-4L s /?o] + 



CH(r« 1 ) 2 + 3C j47c 1 , 



cusp/ 



+ 2r° 



r° 4- ~r 2 r-v ^ 2 

jlhot -t A ^FVTcuspJ 



x[(L h + 2L s )(L h + UL s 



9ir 2 



h,12 L h,21 

A)(7° - 47°) - 4 7 °7° + (7°) 2 + 4(7 S °) 2 , (5.15) 



Bi 



12C\ 



^ f /cusp 



3C 2 A [4L h (c h u - 2L s (j° h - 7j ° + 7s )) + 8L s c h u + 4c[(L h + 2L S 



+ 3 7 ° + 2 7 °) + 4 7 °L 2 - 3^ 2 7 ° - 24 7s °L 2 + 6. 2 7 °)] 
+l2c s u (L h + 2L S ) + Ci/3 (5vr 2 - 6(L| + 4L 2 ))] 



+C F 



r° 

^(24 - r° i21 (L ft + 2L S )) + 7c US p(-^/i + 2L S 



+ ^(2^, 21 c s 12 + ^(7° - 2( 7s ° + A)))) + 2r° M2 4 + (2L S 



Lhfth A2^h ,21 



-4(7° -27° 



ilF f o ^2 
4 v '' 



i cusp / 



[{L h + 2L S ){L\ - 4L 2 + 3vr 2 ) - 36( 3 ] 



+4(7° 



27° 



Aj) - 7°7,% + HlsLh ~ 2 7 °7 S % - 2 7 ° 7 °L S - 7 ?A)^ 



+( 7 °) 2 ^ + 7 ) + 4 7s /3oi> s + 4(7°) 2 i, - 27, 1 , 



(5.16) 



where L 



i, 



ln- 



and L s = In ■ 



2 V~s 



-mF)[i' Explicit expressions for c h , c s and c? are given in the 
appendix. Several comments on this result are in order. Our NNLO expansion are accurate 
to NNLL accuracy, i.e., the coefficients ^-2,1,0 and .64^2,1 are accurate. Complete expression 
for Bq can only be known from a 2-loop calculation, therefore we retain from giving a partial 



result for it here. Similar expansion has been performed in the Ref. [22], but only partial 
NNLL logarithms are presented. In particular, the coefficient of -D 2 and D\ at NNLO is not 
complete in f22]], because of the omission of virtual corrections in their work. We also note 



that due to different definition of S4, in order to comparison with the result in the Ref. [22], 
we have to make the follwing change to j s 



Is 



-Cf7cus P In 



m 



2\2. 



utrrit 



-c F% 



cusp 



hi- 



ts 



m 2 ) 2 \/I 



(u — m 2 )(t — m 2 )mt 



(5.17) 
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We have checked that our expression for #4,3 agree with the Ref. [p2j | after this replacement. 
We also checked that our 11 dependence in the expansion coefficients agree with |2^] for ^2,1 
and -64,3,2 after this replacement. There is an extra \x dependence term in B4 of our expansion, 
CFJcuspLhCii, which is not presented in p^j . There is also 11 dependence in Bq. However 
we are not able to check this term against |22j since the explicit expression for this term is 
not presented there. We give a numerical comparison on the difference of different threshold 
variable definition and the effects of omitting virtual corrections in the resummation result 
in next section. 



6. Numerical discussion 

In this section, we present the numerical results for the threshold resummation effects on the 
s-channel single top production at both the Tevatron and the LHC. The input parameters 
used throughout this section are given below: 



m t = 173.2 GeV, M w = 80.4 GeV, Gp = 1.16639 x 10 

/ 0.9751 0.2210 \ 



-5 



GeV, 



VcKM 



0.2215 
0.0035 



0.9743 
0.0410 



We use the MSTW2008NNLO PDFs M throughout our numerical calculation. 




0.5 1 2 0.5 1 2 0.5 1 2 



Figure 3: Factorization scale dependence of NNLO expanded cross sections for s-channel single top 
production at the Tevatron (1.96 TeV) and LHC (7 TeV and 14 TeV). 

Before presenting the numerical results for resummed cross section, it's important to 
examine to what extent the singular terms approximate the fixed order calculation. It would 
be meaningless if the neglected subleading terms are as important as the singular terms. We 
present the numerical results for the approximate NLO and NNLO cross section, Eq. ( |5.6| ), 
for both Tevatron and LHC in Fig. |3| The factorization scale is chosen as lif = 200 GeV, 
as will be explained below. It can be seen from Fig. ^ that at the Tevatron, which has 
a lower collision energy, the approximate NLO cross section over estimates the exact NLO 
cross section by about 5%. Since the NLO corrections for s-channel single top production 
is quite large, about 40% at the Tevatron, we consider the threshold expansion as a good 
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approximation. Furthermore, the fact that the scale dependence of the approximate NLO 
result is similar to exact NLO result implies the small scale dependence of the subleading 
terms. On the other hand, the NLO approximation doesn't work well at the LHC with 
higher collision energy at 7 TeV or 14 TeV, as shown in Fig. |3[ The differences mentioned 
above are smaller at lower factorization scale and more significant at larger factorization scale. 
Moreover, the scale dependence of the approximate NLO cross section behave quite differently 
from the exact NLO results, indicating that the subleading terms are not only numerically 
large, but also have large impact on the factorization scale dependence of the cross section. 
Therefore, we conclude that our resummation results at the LHC are not as reliable as at the 
Tevatron. Nevertheless, we still give the NNLO approximate and resummed results for both 
the Tevatron and the LHC as a reference. 

In order to calculate the resummed cross section, Eq. ( [4.39| ), we need to determine the 
appropriate scales for the process. In the SCET approach to resummation, there are four 
scales: the hard scale fx^, the jet scale fij, the soft scale /i s , and the factorization scale 
fiF, and we have assumed that the renormalization scale equals the factorization scale for 
simplicity. This is different from the fixed order calculation, where only the factorization 
scale is accessible. This is actually the merit of effective theory, since the calculation has 
been factorized into a series of single-scale problems, and large logarithms can be avoided if 
appropriate value of scale is chosen for each problem. 

First, we choose a default value for the hard scale. Since s-channel single top production is 
similar to the Drell-Yan process, one would expect that the appropriate value of the hard scale 
should be around \f§. However, choosing fj,^ ~ \fi is inconvenient because \fl is a dynamical 
variable. To avoid this inconvenience we set the hard scale at a fixed value, \ih = 200 GeV, 
which is slightly larger than mt. We have checked that invariant mass distribution of the 
virtual W boson peaks around this value, and thus it can be considered as the "average" 
value of \^§. For simplicity, we also choose the factorization scale to be 200 GeV. 

Next, we determine the appropriate soft scale and jet scale. In order to obtain a reasonable 
physical result, we expect that no large logarithms should arise in the soft and jet function 
when the appropriate scales are chosen. In practice, for the determination of soft scale (jet 
scale), we fix the factorization scale at 200 GeV in the factorized cross section Eq. ( |4.39| ). 
and set all the other scales equal to [i s (/ij), and then vary fi s In Fig. ||, we plot the 

cross section which only include the one-loop soft corrections or jet corrections, respectively, 
divided by the tree-level cross section with all scales equal to 200 GeV. From Fig. f|, we find 
that the net corrections are small around fij = 50 GeV for the jet corrections. Thus we take 
it as our default scale choice for jet scale. The choice of soft scale is not as clear as jet scale 
from Fig. |I| since the total soft corrections do not have a clear minimum or maximum. The 
reason is that at NLO, the soft corrections consist of a Drell-Yan like corrections in the initial 
state and a soft gluon corrections in the final state, as depicted in graphs (a), (b) and (c) of 
Fig. |l|. It turns out that the corrections is positive in the initial state and negative in the 
final state. Either initial state corrections or final state corrections show an appropriate scale 
around 25 GeV, but not the sum. On the other hand, an appropriate choice for soft scale is 
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dictated by the picture of underlying factorization, fj, s ~ From our choice for hard and 

jet scale, this implies that soft scale should be chosen around 10 — 30 GeV. In our numerical 
calculation, we have chosen fj, s as 25 GeV, to avoid too close to Aqcd while probe as much 
soft activity as possible. 




ft, (GeV) (GeV) 



Figure 4: Determination of /i s and /ij. On the left (right) is the cross section with only soft (jet) 
corrections, divided by the tree- level cross section, for the variation of /i s (/ij). The red dashed line 
includes only the initial state soft corrections, while the blue dotted line includes only the final state 
soft corrections. The factorization scale is fixed at 200 GeV. 



We also note that the factorized cross section Eq. ( 4.39| ) is derived in the threshold limit, 



S4 —> 0. To capture the non-leading terms, we must match the resummed cross section onto 
the NLO cross section, which can be found in the Ref. Q. We have redone the calculation and 
found complete agreement with the Ref. |J. After matching in momentum space approach, 
the resummed total cross sections is given by 



_RES thres thres 

a = a —a 



NLO 

T <J 



+ (T I . (6.1) 



In table [l], we present the resummed total cross section for s-channel single top and anti- 
top production, as well as the NNLO approximation at both the Tevatron and the LHC. 
All the scales are set to the default values, i.e., fip = = 200 GeV, fi s = 25 GeV and 
fj,j = 50 GeV. It can be seen that the NLO QCD corrections significantly enhance the total 
cross section at both the Tevatron and the LHC ||. The threshold resummation effects 
further increase the NLO cross section by about 3% — 5% at the Tevatron. However, we do 
not observe the large enhancement due to threshold resummation as reported in Ref. |22|| . 
The discrepancy has two origin. First, the 1-loop matching coefficients of hard function was 
not taken into account in |^2|]. second, the definiton of S4 used in does not coincide 
with ours. As explained in Sec. ^, the definition we use also include the effects of collinear 
splitting of final state b-quark, and should be considered as a better choice. We quantify the 
numerical significance of the difference from these different treatment in the end of this section 
briefly. We also show the resummed cross section for single top production at the Tevatron for 
different top quark mass in Fig. |5[ It can be seen that the LO prediction significantly under 
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estimates the total cross section. It's also clear that the resummed cross section dramatically 
improves the scale dependence, comparing with the NLO results. 







0"NLO 


"expand 


CTRES 


Tevatron (top) 


CSlSlffi Pb 


0.443lffl P b 


o.463t8:88i pb 


0.46718:818 pb 


Tevatron (anti-top) 


o.3i8iS:°s P b 


0.4431SS* P b 


o.463i8:88^ pb 


0.46718:818 Pb 


LHC (7 TeV, top) 


2.03l°$ P b 


2.711-- pb 


2.8218:8? Pb 


2.8H8:18 pb 


LHC (7 TeV, anti-top) 


i-i4±8:8i pb 


i-53±8:S pb 


i.6oi8:8I pb 


i-6ol8:85 pb 


LHC (14 TeV, top) 


5.2i+8:i| Pb 


6-91+8:81 pb 


7.i7i8:l Pb 


7.nl8:I pb 


LHC (14 TeV, anti-top) 


3.36™ Pb 


4.46+^ Pb 


4.6418:1° Pb 


4.61183 Pb 



Table 1: Total cross section for single top and anti-top production at the Tevatron and LHC. All the 
scales are chosen at the default values. For the resummed results, The total uncertainties are obtained 
by adding the individual scale variations of (If, A*/i> (J*j, (J-s in quadrature. 




172 173 

m t (GeV) 



Figure 5: Top quark mass dependence of the fixed order and resummed cross section. The band 
corresponds to the variation of factorization scale around 200 GeV by a factor of 2. 



Finally, we give a brief numerical comparision with results presented in Ref. p2[ . As 
mentioned in the last section, our NNLO singular expansion differ from those presented 
in [|22| in two aspects: (a) We have included all the 1-loop matching coeffcients c h ,c s ,c^ in 
our calculation. Therefore our results contain all the NNLL logarithms. 1 (b) The results in 
the Ref. [22] used a different definition of S4, as was mentioned in Sec. [3|. To quantify the effects 
of these differences, we plot two sets of cross sections below. In Fig. [| we plot the NNLO 



1 We refer to Ref. [M] for the accurate definition of logarithmic order in SCET approach to resummation. 
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approximate result with set to zero. As can be seen from the figure, the NLO or NNLO 
approximation fail to approximate the exact NLO result when the virtual corrections are off, 
which give large positive contribution at the NLO. In Fig. ^, we plot the NNLO approximate 



results with the same definition of 54 as in Ref. [22]. The NLO or NNLO approximate cross 
sections are significantly enhanced with such choice for 54, and overestimate the exact NLO 
results by a large amount. 



0.70,- 
0.65- 
0.60- 

2 0.55 

— 0.50 
^0.45" 
b 0.40 
0.35 ■-. 
0.30- 
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NNLO approx(w/o virt.) 
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NNLO approx.(w/o virt.) 

LHC 7 TeV 



— LO 
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— NNLO approx(w/o virt.) 
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Figure 6: Illustration of the importance of the NLO virtual corrections in the NNLO expansion. We 



have set to zero in these plots. 
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Figure 7: Illustration of the effects of different threshold variable definition. We have chosen the 
same 54 definition as Ref. 122] in these plots. 



7. Conclusions 



We have studied the production of s-channel single top quark in the SM at both the Tevatron 
and the LHC. Using SCET, we show that the production cross section can be factorized into a 
convolution of hard function, soft function and jet function in the threshold limit. Each func- 
tion, being sensitive to a single scale, is free of large logarithms once an appropriate scale is 
chosen. By this way, the threshold resummation is performed with the conventional RG equa- 
tion. As a by-product, we obtain a NNLO expansion of threshold singular distributions. We 
also perform a numerical investigation of our resummed formula, using the momentum space 



resummation formalism [3C]. We find that in general, the higher order threshold logarithms 
enhance the NLO cross sections by about 3% — 5% at the Tevatron, and the resummation 
effects significantly reduce the factorization scale dependence of the total cross section at 
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the Tevatron, while at the LHC the factorization scale dependence has not been improved, 
compared with the NLO results. 
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A. Relevant anomalous dimensions and matching coefficients 



The various anomalous dimensions needed in our resummation can be found, e.g., in the 
Refs [|l|, H 0- We list them below for the 

convenience of the reader. The QCD j3 function 

is 

2 



with expansion coefficients 



-2a, 



_ a s _ / a s \ A 



(A.l) 



h = yCl - Y° ATFn f ~ 4C F T F n f , 



2857 



C\ + I 2C F —CfCa 



205 



1415 



Ci) T F n f + {^C F + ^§C A ) 2fnJ,(A.2) 



where Ca = 3, Tp = 1/2 for QCD, and nt is the number of active quark flavor. 
The cusp anomalous dimension is 



7cusp(a s ) = 7cus P J + Tcusp (f^) + • • • , 



(A.3) 



with 



7° 
/ cusp 


= 4, 




^cusp 


= 4 


Y67 


7 2 
/ cusp 


= 4 


~c\{ 



3 y Ca " Y TFnf 



27^ + 45^ + ~ Cs ) + 



22 
T 



418 40 2 56 

yr + 27 n ~ y^ 3 



+C>2>n/ ("f + 16C 3 )-^nJ 



(A.4) 
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The other anomalous dimensions are expanded as Eq. ( A. 3 ), and their expansion coeffi- 



cients 


are 


7,° = 


— 3Cp, 


l\ = 


c%{- 


ll = 


-2C F , 


7q = 


CfCa 


7^ = 


3C F , 


7^ = 


<*(§ 


7? = 


— 3Cf, 


7) = 





+ 2tt 2 - 24C 3 + C F C A 



961 11 2 



54 6 



^ + 26C 3 ) + CVTpn/ ( — + -vr 



130 2 



98 _ \ 40 ^ m 



17 22 



2^ + 24C 3 ) + C F C A ( — + — - 12C 3 ) - CTpTjrn/ ( - + ) , 



2 8 



3 9 



+C F T F rif 



+ 2tt 2 -24C 3 1 +(>(" v 
'242 4 



1769 11 2 , 

vr 2 + 400= 

54 9 V 



27 9 



-7T 



(A.5) 



7h and 7 S can be obtained from the anomalous dimensions above through the following 
equations: 



Ih = ?>lq +7Q, 

7 S = -274,-7^ + 7^. (A. 6) 

The hard function is a 2 x 2 matrix in color space. To 0(a s ), it can be written as 



H 



1 _|_ 2<L H l ^LfJ 1 



as. ttI 
47r jfl 21 







(A.7) 



can be obtained from evaluating the first two diagrams of Fig. [2] and the corresponding 
counter-terms. It is given by 



(A.8) 



with 



-11 



C F 



2 . - 9 a t In — — + ln(l - x t )(-2x t + 8 In x t + 6) - 6 ln 2 (l - x t ) 
t — mi 



-4Li 2 



at - 1 



1 - x t 

29 

4 In 2 x t + (2xt - 12) In x t + — vr 2 - 28 

6 



(A.9) 



where xt = m\js. We have checked this result against the existing NLO virtual corrections to 
s-channel single top production || and found complete agreement. The remaining diagrams 
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of Fig. H corresponds to H\ 2 and H\\- They do not contribute to the cross section at the 
NLO but do at the NNLO, therefore are required for a complete NNLL resummation. The 
results are 

with 

= -2 [m 2 (s + 2t) + M^s - s 2 - 2si - 2p] C (0, 0, a, M^, 0, 0) 

+J^l [ m * (*( 3 * - 2M w) + * 2 + ^i) + (s + i) 2 (M^ - a - 2t)] C (0, u, m 2 , M^, 0, 0) 

+4i(m? - t)C (0, 0, s, 0, 0, Myy) - u{M^ + u- i)C (0, 0, u, 0, 0, 0) 

+{m 2 - s)(M^ + u- i)C (0, S, m 2 , 0, 0, M^) - u(M^ +u- 4)0,(0, u, 0, M^, 0, mf ) 

+(s + t)(M& + u - i)C (m 2 , 0, «, m 2 , 0, 0) + 2m?*Co(0, t, 0, M^, 0, m 2 t ) 

} ; [mf + m^M^ - 3s - 2t) + m 2 (3s(s + 2f) - 2M$ v {s + t)) + s 2 (M^ -a- 2t)] 
mf — s 

xC (m 2 , S, 0, m 2 , 0, Af&r) - 2t 2 (m 2 - i)D (0, 0, 0, m 2 t , i, a, 0, 0, 0, Afjy) 
2t(m 2 - £) 2 A)(0, s, 0, t, m 2 , 0, m 2 , M^, 0, 0) 

+(s + t) [m 2 (M^ -a- 2t) + - 2M$ v (s + t) + s 2 + 2si + 2t 2 ] 
xD (m 2 , 0, 0, 0, «, a, m 2 , 0, 0, M^) 

- mf(M^ -a- 2i) + m 2 (M^ - 3M^(s + £) + 2s 2 + 5si + At 2 ) 

-M^(s + t) + 2M$ v (s + £) 2 - s 3 - 3s 2 £ - Ast 2 - 2t 3 D (0, 0, 0, m 2 , s, u, M$y, 0, 0, 0) 

+ 7 5 * A (£o(m 2 ,0, M 2 ,) + 5 (m 2 , 0, m 2 )) + -^-^B (s, 0, M 2 ,) 
(m/ — a)(a + t) mf - s 



^iu /* s „ /* 9,> 2 tit — mf), u(u — mf) 

-JB (u, 0, 0) + B (tt, 0, m 2 )) . V : In (A.ll) 



a + t K u ^ ' ' ' ' u ^ ' ' *" e m s-M 2 , i(i-m 2 ) 

where the Bq, Cq and Dq are the conventional Passarino-Veltman function |70|1 , evaluated at 
the point where the 't Hooft mass is set to mt- For example, the Bq function reads 

B^M) - J {p _ ml + m \ pi? _ ~, + , £j ■ (A.12) 

and similarly for Co and Do. This is just a simple way to extract the /j dependece from -ffj^- 
The analytical form of the singular Passarino-Veltman function can be found, e.g., in the 
Ref. J7lj. Note that H\ 2 and H ^ are UV and IR finite. Given the 1-loop matching coefficient 
above, one can check that the RG evolution equation of the hard function exactly has the 
form of Eq. flPop . 

The calculation of the soft function can be divided into the calculations of soft integral Is 
and the corresponding color factor. We first discuss the soft integral below. The soft integral 
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corresponding to diagram (a) of Fig. |] and its mirror image counterpart can be written as 



la ~ ^ 1 4vT 



n a ■ n b 



<?)(«& • q)' 



(A.13) 



where we work in n = 4 — 2e dimension, and the factor of 2 comes from doubling the contri- 
bution of diagram (a) by including its mirror image counterpart. This integral is evaluated 
by Becher and Schwartz j47p, with the result: 



la 



4tt 



2 In 



2 2n ab _ if_ 
n la n lb 3 



5(k + ) + 16 



1 

k+ 



In 



/ 2n a • n b 



1 1 



[*+■**]' 



(A.14) 



where n 



a(b) 



n 



a m • n\ and [/]i a,b ' is the star distribution defined in the Ref. (45]. Note that 
we have put a bar on / to denote that divergent terms have been subtracted in MS scheme. 
The soft integral corresponding to diagram (b) and (c) of Fig. || reads 



h 



2 ( 



47T 



h = 2g. 



Ait 



■S(q 2 )e(q )5(k + -n l -q)^—^. 

/d n n 



(ni • q)(v ■ q)' 



(A.15) 
(A.16) 



where v + = v ■ n\ . The simplest way to do the integral of I b and I c is working in the lightlike 
coordinates along the n\ direction, in which any four vector can be written as 



a „ — U i + -M i A 1 



The results for the integrals in Eqs. ( A.15| ) and ( [A. 16 ) are 



(A.17) 



I b = _i { 41nv + (5(A; + ) -4 

47T 



[fc+, M ] ' 



47T 



-4 In 2 tf 



7T 

~6~ 



6{k + ) 



1 In 



[fc+, M ] • 



(A.18) 
(A.19) 



The next two diagrams, (d) and (e), vanish, as was explained in the Ref. [47]. Diagrams 
(f) and (g) are the most complicated diagrams to be evaluated. They read 



if = ig. 



2 ( A 7g 

47T 



j^-^Zj6(q 2 )0(qo)5(k + - m ■ q)- -A— 

(2ir) n 1 (n a ■ q)(v ■ q) 



4ir 



d n q 



n b ■ v 



77y )6(q )5(k + - m ■ q)- v -. 

(27r) n 1 (n b -q){v-q) 



(A.20) 



(A.21) 



We calculate diagram (f) first. As before, we work in a lightlike coordinates along the n\ 
direction. The delta function in Eq. ( [A. 20 ) can be used to integrate out dq + and dq±. After 
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some simplification we arrive at 



a s (pv+) 2 * 4e^ 



4tt (fc+) 1 + 2e 



Vir(i/2 - e) 



dx 



(19 



x 



;,, sm 2t (1 + z) (x + p a — 2^Jxp a cos 6) ' 

(A.22) 

where we have defined p a = — \ n " , a; = 9 fe: ^ , and 6 is the angle between q± and n a j_. 

Note that we have (1 + p a ) = 2n a ■ v^. A somewhat similar integral appear in the evaulation 
of hadronic thrust distribution [fT2j] . The integral in Eq. (A.22) is straightfoward to do, with 
a result 



\2c 



4tt (k 



+ U+2e 



-- + 81n(l + p a ) + e(8Li 2 ( 

e V V 1 + Pa 



7T 



2\ 1 



41^(1 + p a )- — 



(A.23) 



Using the expansion of {k 



+ \-l-2e 



in terms of star distribution [45]: 



J_ (JL\ 2 e 
k+ \k+) 



l 

2e 



k+ 



[k+,ti] 



- 2e 



1 , k+ 

7T m 



+ 0(e 2 ), 



(A.24) 



we obtain 



If — — 



a. 
-in 



+8 



-4Li? 



Pa 



1+Pa 

fc + 1 + Pa 
fJL V + 



+ 21n 2 (l + p a ) 



8 In v + ln(l + p a ) + 4 In 2 v + H 



<5(£; H 



(A.25) 



Diagram (g) can be obtained from diagram (f) by the simple replacement I g = —If(p a — > p^). 
Given the integral above, the soft function can be derived by combining them with the 
corresponding color factor: 



T—Sll = C\CFiJa + h + h), 
47T 

^^12 = \ C AC F {Jf + Ig), 



D 21 — D 12; 



7 21 — "12 
4^ 22 " 4 



— ^22 = —-;CF(I a + h) + -CAC F h + -Cf(CX — 2)1 f — -Gpl g . 



(A.26) 



We note that the complete set of color factor for all 2 — > 2 processes has been worked out 



in the Ref. [72]. To check that the soft function we obtain indeed obeys the RG evolution 
Eq. ( [4.34 ), it's convinient to make a Laplace transformation to the soft function. Explicitly, 
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we have 



s l n {L,n) = C\C F ( 4L 2 + 8 In 2Ua ' "f + L-4L) .- . 



l + Pb 



ntn^ 



'12; 



s\i{L,p) = s\ 2 (L,fi), 

sl 2 (L, M ) = C F (C\ - 1)L 2 - 2C A C F L - 2C F In 

+2C F {C 2 A - 2) In i±^L + 4C F In + c| 2 , 



(A.27) 



where L = In - . and 



7T 



= C7iC F (21n 2 ^^-41n 2 U + +41n V + + 

V "la^ib 6 



12 



CaCf —2Li 2 



Pa 



+ 2Li 2 



P& 



.i + Po/ Vi + Pfc 

+ ln 2 (l + p a )-ln 2 (l + p 6 )Y 



c! 2 = CfCI -Li 



41nt>+ ln(l+p a ) +41nu + ln(l + p 6 ) 



21n^n(l + p a )+ln 2 «+ + Jln 2 (l+p a ) + ^7r 2 
1 + p a / 2 24 



Pa 



+c F 



1 



In 2 



2ra ab 
niani h 



+ 2Li 2 



Pa 



1 + Pa 



2Li 2 



1 + Pfe 



+ 41nv+ ln(l + p a ) 



7T- 



-4 In v + ln(l + p b ) + ln z v + - ln"(l + p a ) + ln z (l + p 6 ) - —j + 2C A C 2 F lnt; + (A.28) 

Using the expressions above, we confirm that the RG equation of the soft function agrees 
with Eq. ( [4,34 ) . This shows that our resummed cross section is RG invariant in the threshold 
limit, which can be considered as a non-trivial check of our result. 
Finally, the Laplace transformed jet function is given by [frSfl: 



j(L,p) = l + 



«s(p) 

4-7T 



-C Fl ° cusp L 2 +^L + c ! 1 



(A.29) 



with cj = C F x (7- §tt 2 ). 
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